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ABSTRACT 

We discuss the constraints that future photometric and spectroscopic redshift 
surveys can put on dark energy through the baryon osciUations of the power 
spectrum. We model the dark energy either with a perfect fluid or a scalar 
field and take into account the information contained in the linear growth 
function. We show that the growth function helps to break the degeneracy in 
the dark energy parameters and reduce the errors on wq , wi roughly by 30% 
making more appealing multicolor surveys based on photometric redshifts. 
We find that a 200 deg^ spectroscopic survey reaching z « 3 can constrain 
wq.wi to within Awo = 0.21, Awi = 0.26 and to Awq — 0.39, Awi = 0.54 
using photometric redshifts with absolute uncertainty of 0.02. In the scalar 
field case we show that the slope n of the inverse power-law potential for dark 
energy can be constrained to An = 0.26 (spectroscopic redshifts) or An = 0.40 
(photometric redshifts), i.e. better than with future ground-based supernovae 
surveys or CMB data. 

Key words: 

1 INTRODUCTION 

Dark energy is the most elusive component of the universe: it is dark, it does not cluster and, 
in most models, it is subdominant in cosmic history until recently. On the other hand, it is 
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also the most pervasive component, in that it affects baclcground expansion, hnear dynamics 
and possibly non-linear as well. 

The characterization of DE has been so far based almost uniquely on background tests 
at rather low redshifts {z < 1: Riess et al. 1998, Perlmutter et al. 1999, Tonry et al. 2003, 
Riess et al. 2004) or very large redshifts {z ~ 1000: e.g. Netterfield et al. 2002, Halverson 
et al. 2002, Lee et al. 2002, Bennett et al. 2003), although tests at intermediate redshifts 
based on strong lensing have also appeared (Soucail, Kneib & Golse 2004). These tests are 
based essentially on estimations of luminosity Dl or angular-diameter distances D{z), i.e. 
on integrals of the Hubble function H{z) which, in turn, contains an integral of the equation 
of state . If we aim at making real progress it is then necessary to cross this information with 
probes at intermediate z and with other observables such as the linear growth function G{z). 
The goal of this paper is to show how future observations of large scale structure (LSS) at z 
up to 3 can set interesting limits on models of DE taking advantage of the complementarity 
between H, D and G. 

The method we use is based on recent proposals (hinder 2003, Blake & Glazebrook 2003, 
Seo & Eisenstein 2003) to exploit the baryon oscillations in the power spectrum as a standard 
ruler calibrated through CMB acoustic peaks. In particular, Blake & Glazebrook (2003) and 
Seo & Eisenstein (2003; hereinafter SE) have shown the feasibility of large (100 to 1000 square 
degrees) spectroscopic surveys aX z 1 and z ^ ?> to put stringent limits to the equation 
of state w{z) and its derivative. SE have also shown that the redshift error introduced by a 
photometric survey (with a relative error on 1 -|- 2; of 1%) could be compensated only by a 
20-fold increase of the sampled volume. 

This paper has a threefold motivation. First, we extend the Fisher matrix method by 
taking into account the information on the DE evolution contained in the growth factor G 
at various redshifts. While SE marginalized over G, we find that its marked dependence 
on the cosmo logical parameters, especially on fi^, associated with its different direction of 
degeneracy, are helpful in narrowing down the constraints. Secondly, we extend the analysis 
to a scalar field model of DE with inverse power-law potential (Ratra & Peebles 1988), based 
on the inverse power-law potential: this allows us to evaluate the growth function without 
further assumptions and to exploit the properties of tracking solutions (Steinhardt et al. 
1999). Third, we pay special attention to the possibility of testing DE with a photometric 
redshift survey. Indeed the effect of the growth factor in reducing the confidence regions of 
the cosmological parameters makes more appealing the use of photometric redshift surveys 
to show any possible redshift evolution of the DE. We explore the constraints derived on the 
cosmological parameters from relatively wide multicolor surveys where the absolute redshift 
uncertainty keeps almost constant around 5z = 0.02 — 0.04 in the redshift range z = — 3.5. 
Such surveys can be realized with the current wide field imagers on 4-8m class telescopes. 

Finally, it is interesting to speculate on another aspect of the method. The baryon oscil- 
lation scale measures the angular-diameter distance D{z) while the SNIa method observes 
the luminosity distance. As remarked recently by Bassett & Kunz (2003) the two quantities 
are related in a trivial way only in standard cosmologies: new physics like photon decay can 
lead to an observable breaking of their "reciprocity" relation. A photometric survey can be 
adapted to the scope since the degradation of the error in D{z) due to the redshift smearing 
is relatively reduced (with respect to the error in H{z)), as it will shown below. 
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2 FISHER MATRIX METHOD 

The method proposed in SE is based on the Fisher matrix, an approximation to the hkehhood 
function that provides under some conditions the minimal errors that a given experiment 
may attain (Eisenstein, Hu & Tegmark 1999, EHT ). Here we review the method adopting 
the notation of SE, to which we refer for further details. Our starting point is a cosmological 
model that predicts the evolution of the Hubble parameter, H{z), and the growth of linear 
perturbations, G{z). From H{z) we can estimate the angular-diameter distance in flat space 
as ( we put c = 1) 

Let us choose a reference cosmology, denoted by a subscript r, for instance ACDM, and 
calculate the matter power spectrum in real space at z = 0, Por(fcr)- Then the predicted 
observable galaxy power spectrum in a different cosmology C (no subscript) at z in redshift 
space is 

P^UZ; K, t^r) = Ps{z) + + Pl^rPorik) , (2) 

where k is the wavenumber modulus and yU its direction cosine. Several comments will clarify 
the meaning of this equation. Ps{z) is a scale-independent offset due to incomplete removal 
of shot-noise. The factor {D^ / H) / {D"^ / Hr) takes into account the difference in comoving 
volume between the two cosmologies. The factor (1 + (3^^Y is the redshift distortion. G{z) 
is the growth factor of the linear density contrast 5mi 

G{z) = 5„,{z)l5o, (3) 

where 5q is the present density contrast (and therefore G{z) is normalized to unity today). 
The density parameter VLm{z) depends on the cosmology. In fiat-space ACDM it is given by 

= „ , tf o ■ w 

ilmo[^ + ZY + 1- ilmO 

and it is therefore parametrized by the present density parameter Vtmo- The bias parameter 
I3{z) = Q'!^{z)/b{z) (assumed scale-independent) is evaluated for the reference cosmology 
using the formula 



^s,, = m-'^t%z)as,Uz)^l + ^ + m , (5) 

where ag is the variance in spherical cells of 8 Mpc/h for galaxies (subscript g) or matter 
(m) (notice that the variance as is calculated using the power spectrum averaged over fi). 
Finally, we can include a redshift error by rescaling the power spectrum 

P = Potsc-'"^"''" , (6) 

where 

ar = {5z)/H{z) , (7) 
is the absolute error in distance and 5z the absolute error in redshift. 
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We need also the relation between the line-of-sight and transverse wavenumbers k\\ 
and the reference values: 

= k^D/Dr, (8) 

kr\\ = k\\Hr/H. (9) 

From these relations we derive the relation between the wavenumber modulus k and the 
direction cosine /i = k ■ r/Zc in the reference cosmology and in the cosmology C 

k = Rkr, (10) 

^ = TTr' ^^^^ 

where 



^= DHr • ^^^^ 

The observed power spectrum in a given redshift bin zq depends therefore on a number of 
parameters, denoted collectively pt, such as flmo, G{zo), H{zq) etc, as detailed in Table I. The 
redshift-dependent parameters [Pg, H,G, D, (3 ) are assumed to be approximately constant 
in each redshift bin. Then we calculate, numerically or analytically, the derivatives 

Finally, the Fisher matrix is (Tegmark 1997) 

where 

nP{k, fi) 



Veff{k,fl) 



^survey ; (-^^) 



nP{k,n) + l_ 

n being the number density of galaxies. 

In the following we will use also the CMB Fisher matrix (Seljak 1996, EHT) 

Fij,cMB = 1^1^ ~fi;r^^'''''^^y^hr ' ^^^^ 

e x,Y '-^t^j 

where Cx/ is the multipole spectrum of the component X, and X = T,E,B,C denote 
the temperature T, E and B polarization, and the cross-correlation TE, respectively, with 
covariance matrix {Cove)xY- We adopt the same specifications used in EHT and in SE, cor- 
responding to an experiment similar to the Planck satellite. The spectra have been obtained 
with the Boltzmann code Cmbfast (Seljak & Zaldarriaga 1996). Once we have the Fisher 
matrices for all data we sum to obtain the total combined matrix 

-^ij Pij,CMB ~l~ ^ ^ Pij, survey ■ ^^^) 

The independent parameters to be included in the Fisher matrix are listed in Table I. 
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Table I. Parameters 


Poik) 


1 


reduced total matter Um 




2 


reduced baryon dens. Ub 




3 


optical thickness 


T 


4 


primord. fiuct. slope 


ris 


5 


total matter dens. 


^mO 


For each survey at z 


6 


shot noise 


Ps 


7 


ang-diam. distance 


\ogD 


8 


Hubble param. 


log if 


9 


growth function 


logG 


10 


bias 


log/3 


additional CMB parameters 


11 


decoupling ang.-diam. dist. 


\ogD 


12 


tensor-to-scalar ratio 


T/S 


13 


Ci normalization 


log A, 



The main difference with respect to the method of SE is that we exploit the information 
contained in the growth function. As it will be seen below, this information is of great help 
in reducing the bounds on the DE parameters, for two concurrent reasons: first, G has a 
marked dependence upon the cosmological parameters, i.e. dlogG/dflmo is at 2; > 1 as large 
as dlogD/dflmo and quite larger than dlogH/dQrno (see Fig. 1); second, the parameter 
degeneracy in G{z) has a slightly different direction with respect to the degeneracies in H 
and D (see the discussion in Cooray, Huterer & Baumann 2004 ). 
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Figure 1. Derivatives dlo^X/dVlmO as function of 2, where X = H,D,G. 



Unless differently specified, the reference cosmology is AC DM with the same parameters 
as in SE, = 0.35, h = 0.65, Uf, = 0.021, r = 0.05, = 1, T/S = 0, = . 



3 THE SURVEYS 

One of the main aim of the paper is to show how it is possible to build, using deep and 
wide multicolor surveys, a photometric redshift sample of galaxies suitable for the analysis 
of the galaxy power spectrum as a function of redshift. Previous analyses reached two main 
conclusions about the possible strategy. First of all spectroscopic surveys are clearly favorite 
with respect to photometric surveys because any uncertainty > 1% in the photometric 
redshifts increases the uncertainty on H{z) a.t z > 0.5 to > 10%, for surveys covering 
areas of 100 — 1000 deg^. Second, the relatively bright magnitude limit needed to keep a high 
completeness level in the spectroscopic follow up implies a low number density of sources 
per unit area. Large survey areas are hence needed to guarantee adequate total number 
of objects. Here we show that it is possible to envisage a photometric redshift survey with 
adequate depth and limited area of the sky at 2; = 0.5 — 1 and at 2; = 2. 7— 3.5 that can provide 
constraints on the redshift dependence of the DE equation of state which are competitive 
with those obtained from SNIa and CMB in the foreseeable future. 
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To measure the redshift evolution of H{z) and of the angular diameter it is important to 
analyse the power spectrum in a wide redshift interval with specific sampling at z = 0.5 — 1 
and especially at z = 3 where the linear regime of galaxy clustering extends to smaller 
physical scales allowing the analysis of acoustic oscillations to higher frequencies in the 
range 0.1 < A; < 0.5 h Mpc~^. What is important is the estimate of the characteristic 
wavelength of the acoustic oscillations, kA = 27r/s which is connected with sound horizon s 
at recombination. For our standard ACDM model kA — 0.06 h Mpc~^. It is important to 
emphasize that local surveys could probe only the first peak at ~ 0.075 h Mpc~^. Thus, 
an appropriate comparison with the CMB power spectrum requires higher frequencies and 
consequently high redshift surveys, in particular at z = 3. 

Indeed, large scale galaxy surveys at 2; > 3 could in principle be built in a straightforward 
way using the standard multicolor dropout technique which exploits the detection in the 
galaxy spectra of the UV Lyman absorption due to the local ISM in the galaxy and in the 
intergalactic medium present along the line of sight. At z = 2.5 — 3.5 the Lyman series 
absorption progressively enters the UV filter causing a reddening in the U-R galaxy color. 
At z = 3.5 — 4.5 the dropout shift in the blue causing a reddening in the B-R galaxy color and 
so on. Spectroscopic confirmation has provided a relatively high success rate 75% (Steidel et 
al. 1996, 1998) although in practice it is difficult to obtain a high completeness level in the 
magnitude interval R = 24 — 25.5, which is the present limit for systematic spectroscopic 
follow up at 10m class telescopes. A survey completeness level of 90% could be reached more 
easily at -R < 24.5 but with a surface density of ~ 0.3 arcmin"^ at z = 2.5 — 3.5. The color 
selection of the sample requires very deep exposures in the UV band implying efficient UV 
imagers at 8-lOm class telescopes to cover a large area in a short time. Surveys at 2; > 3.5 
which use a B-dropout technique can avoid high UV efficiency but the surface density of 
objects at the spectroscopic limit of / ~ 25 is only ~ 0.1 arcmin"^. In other words galaxies 
at z ~ 4 — 5 are rarer and fainter. With these considerations, surveys at z ~ 3 down to a 
magnitude limit of i? ~ 25 are the best compromise for the study of large scale structures at 
high redshift provided that efficient wide field imagers with high UV efficiency are used. The 
Large Binocular Camera in construction at the prime focus of the Large Binocular Telescope 
is an ideal imager to this aim (Pedichini et al. 2003). 

It is well known from the previous studies that there are two sources of statistical error in 
the power spectrum estimate. The sample variance, due to the number of wavelength bins, 
depends on the survey volume and on the scale of the selected acoustic peak in the linear 
regime (e.g. Peacock & West 1992; Blake & Glazebrook 2003): 



^ P ^ V 47rk^Ak ' 

The shot noise is the other source of noise. For a fixed volume it decreases with increasing 
nP > 1, where n is the comoving galaxy number density. In practice it is possible to show 
that the decrease of the shot noise saturates for nP » 1 and values in the range 1 < nP < 5 
are more appropriate (SE). 

To reach at least the second peak at A; = 0.135 h Mpc~^ it has been shown that adopting 
AA; = kA/4: = 0.015 h Mpc~^ to ensure an adequate wavelength sampling, a value ^ = 0.02 
is needed as the fractional height of the peak is about 4% (Blake & Glazebrook 2003). This 
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requires that volumes of the order of 3.6 x 10* Mpc^ should be mapped at z = 0.5 — 1 
and at 2; = 3. Since the volume per square degree is ~ 3.5 x 10^ h~^ Mpc^ deg~^ at 
z = 2.5 — 3.5 and = 1.7 x 10^ Mpc^ deg~^ at z = 0.5 — 1.3 this implies surveys at least of 
the order of 100 — 200 deg^ at 2 ~ 3 and 2; ~ 1, respectively. 



Concerning the shot noise, the condition of 1 < nP < 5 implies n ~ 4 — 20 x 10 '^u'l^K^ 
Mpc"^ respectively, given P k. 2500a|^g (Mpc//i)3 at A; = 0.2 /i/Mpc. 



At photometric redshift survey down to the magnitude limit of i? < 25 gives a 

total number density in the range z = 2.7 — 3.5 of n ~ 2 x 10~^ Mpc~^ after integration 
over the UV Schechter luminosity function of Poli et al. (2001). 



At z = 1, the same photometric survey at i? < 25 in the range 2; = 0.5 — 1 would give a 
total number density of ~ 4.5 x 10~^ Mpc~^ as obtained from the blue luminosity function 
of Poli et al. (2003). For sake of comparison, we adopt the same redshift binning at 2 ~ 1 as 
SE. 



Detailed densities and volumes for the basic surveys of area lOOdeg^ are shown in Table II 
in the same redshift intervals used by SE (except we restrict conservatively the deepest survey 
to z > 2.7 instead of z > 2.5). We will present results also for larger surveys, corresponding 
to areas of 200 deg^ and 1000 deg^ respectively, and with sparser surveys at 2 ~ 1, adopting 
a density 10 times smaller as that quoted in Table II. The absolute error in redshift is taken 
to be either zero (spectroscopic survey) or 6z = 0.02 and 6z = 0.04 (photometric surveys). 
Finally, we also include results for a hybrid survey, i.e. a photometric survey with 6z = 0.02 
at z ~ 1 combined with a spectroscopic at 2; ~ 3. This choice depends on the fact that multi 
object spectroscopy of a complete color selected sample of Lyman break galaxies at 2; = 3 is 
able to give more stringent constraints on the cosmological parameters (see Fig. 3) than a 
similar survey at z ^ 1 in the same area. Moreover a spectroscopic survey at z ^ 3 appears 
also more efficient in terms of multiplexing and completeness level if the magnitude limit 
is i? < 25, corresponding to a surface density < 1 arcmin^. At 2; = 1, on the contrary, the 
average surface density even at the relatively shallow limit i? = 23 is already of the order of 
8 arcmin^ implying a multiplexing of more than 460 targets in a field of e.g. 7x8 arcmin^ to 
be compared with the multiplexing capabilities of recent instruments at 8m class telescopes 
like VIMOS at ESO/VLT where only 150 targets can be observed in a single exposure in the 
same area. 
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Table II. Surveys 




z 


Vs (Gpc3) 




n 


0.5 


-0.7 


0.030 


6.9 


■10-2 


0.7 


-0.9 


0.040 


4.2 


■10-2 


0.9 


- 1.1 


0.050 


3.1 


■ 10-2 


1.1 


- 1.3 


0.057 


2.4 


■ 10-2 


2.7 


-3.5 


0.27 


2 ■ 


10-3 



As in SE, to these surveys we add data from SDSS (Eisenstein et al. 2001). The SDSS 
data is always considered spectroscopic. 

Redshift errors of the order of 1-2% from photometric surveys are demanding but still 
reachable if the photometric accuracy is 6m ^ 0.02 and if a good wavelength coverage and 
sampling is guaranteed (see e.g. Wolf et al. 2001). This requires broad and intermediate band 
imaging on wide field imagers at 8m class telescopes like Suprim-Cam at the SUBARU or the 
upcoming Large Binocular Camera at the Large Binocular Telescope (LBT). In particular 
the use of intermediate band filters in the optical range can be important to ensure the 
required redshift accuracy in the specific redshift intervals 2; ~ 1 and z ~ 3. In the redshift 
interval z = 1.5 — 2.5 near infrared imaging is needed on wide areas but this requires a new 
generation of NIR instruments. 

To the redshift surveys we add CMB data, adopting the same Planck-like specifications 
of SE and EHT. Finally, we evaluate SN constraints from surveys that reproduce approxima- 
tively near future catalogs from ground-based observations: 400 SNIa distributed uniformly 
in redshift between z = and z = 1.5, with magnitude errors Am = 0.25 and the same 
stretch factor correction as in the sample of Perlmutter et al. (1999). We marginalize over 
the i?-dependent magnitude Ai and over the stretch factor. These SN likelihood functions 
are therefore completely independent of the Hubble constant determination. We also com- 
pare with SN likelihood contours obtained assuming an independent 10% gaussian error on 
Hq. Notice that the SN likelihood contours are estimated through the full likelihood function 
rather than with a Fisher matrix to facilitate the comparison with current estimates. We 
will compare the LSS forecasts also with the present SN dataset of Riess et al. (2004). 

For the bias of the galaxy populations at z = 1 and z = 3 we adopt the same values as 
in SE, i.e. 6 ~ 1 — 1.5 at 2; ~ 0.5 — 1 and b ~ 3.3 at 2; ~ 3 (see Adelberger et al. 1998). 
Specifically, we assume as^miz = 0) = 0.9 and crs,m{z) = a8,m{0)G{z) and as,g = 1 at 2 = 3 
and (Tg^g = 1.8 at all other redshift bins. 
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4 PERFECT FLUID DARK ENERGY 

To review notation and to check our code against the results of SE we begin the analysis 
with their DE model, based on a perfect fluid with equation of state p = w{z)p where 

w{z) = Wo + wiz . (18) 

Notice that such a perfect fluid model can be only a phenomenological description of the 
background dynamics at low z. Beside the obvious fact that w diverges at high z, the problem 
is that a meaningful derivation of the fluctuation equations requires as additional prescription 
the dependence of w on the perturbed dark energy density. Moreover, if one assumes for 
simplicity that 6w/6p = 0, then dark energy becomes highly unstable at small scales since 



the sound speed Cg = ySp/6p is imaginary for w{z) < 0. We can discuss therefore the growth 
function only in the ACDM limit Wq = —l,Wi = 0. In this limit a useful approximation is 

dlogG 

— — = nm{ay, (19) 

a log a 

with 7 ^ 0.6. Since G can be calculated only in the ACDM limit, we assume that G does 
not depend on wo,Wi: this is a conservative statement, since any dependence would further 
restrict the area of the confidence region. Moreover, we found that, naively extending the 
validity of (fT^ to other values of Wq,Wi (i.e. assuming dark energy to be completely homo- 
geneous), the results do not change sensitively. Therefore G will be assumed in this section 
to depend only on Qmo- The derivative with respect to flmo is 

dlogG f ..rflogOn f e^-nUc^r 

— — = 7 / i2m(a)^— da = — 7 / — — - — t-- — da. (20) 

Once we have the total Fisher matrix we invert it and extract a submatrix Fsub,ij con- 
taining all the parameters Pi that depend on the cosmological parameters u!m,^mo,Wo,Wi. 
Beside Um, ^mo these are, for each survey, 

logD,\ogH,logG, (21) 

plus log D for the CMB. The diagonal elements of Fg^l ,j^j give the square of the errors on the 
relative variables. The relative errors on H,D,G are plotted in Fig. 2 for the various surveys 
and redshift errors. These intermediate results are consistent with those given in SE (who 
however considered a different combination of surveys and volumes). As remarked in SE, the 
errors on D (contrary to H and G) depend only weakly on redshift uncertainty. At 2; = 1 
the angular diameter distance can be estimated to 4.3 % in a spectroscopic survey of area 
200 deg^ and to 7.5%(10.0%) in a similar survey with 6z = 0.02(0.04). Roughly a factor of 
2 is gained by extending to 1000 deg^. 
Changing to the DE variables 

dn = (^^mO^^ ^mO, U!o, ^l) , (22) 



we obtain finally the projected DE Fisher matrix 
_(dpA ( dp J \ 



(23) 
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Then we again invert and take a submatrix M^jg^wi whose eigenvalues and eigenvectors de- 
fine the confidence elhpsoid on the plane Wq,Wi. Notice that a two-dimensional confidence 
region corresponding to la includes only 39% of likelihood. Here and in all subsequent plots 
the confidence regions are given to 68% : this requires therefore that the la values of the 
parameters (i.e. the eigenvalues of the two-dimensional Fisher submatrix) are multiplied by 
1.51. 

The results for the various surveys and for different redshift errors are in Figs. 3-6. In 
Fig. 3 and 4 we show the effect of marginalizing over G versus fixing G = G{Qmo) as in 
Eq. (fT^ for the four surveys at 2; ~ 1 and the survey at 2; = 3, always including CMB 
and SDSS. As it can be seen, the effect is rather strong, especially at z = 3, where the 
information in G almost halves the errors in wo,wi. At this redshift, we could say that a 
survey of 200 deg^ that includes the information in G (from now on, we call this approach 
"including G") is roughly equivalent in terms of constraining power to a survey of 1000 deg^ 
with G marginalization ("marginalizing over G"). The effect at 2 ~ 1 is smaller, because 
dlogG/dflrno is smaller at low z (see Fig. 1). 

In Fig. 5 we show the results for all redshift bins marginalizing over G as in SE while 
in Fig. 6 we include G{z). In the same plots we show the contours of the likelihood for 
simulated SNIa data. These are obtained by a full likelihood analysis of simulated datasets, 
marginalizing over the Hubble parameter derived from the SNIa themselves. This renders 
SNIa sensitive to the slope dlogD^/dz of the luminosity distance D^, rather than to its 
absolute value. Since 

^^"g-P^ oc - (24) 

dz H{z)Dl{z) ' ^ ^ 

it is clear that the SNIa data is nearly degenerated in the direction of HD = const, at 2; ^ 1. 
In Figs. 5-6 we also show the 68% confidence region for SNIa obtained assuming instead a 
10% gaussian error in Hq. 

The la errors on Um, ^m, wo, wi are reported in Table HI and IV, where we include always 
G{z). For Um we give the relative errors; for the other parameters the absolute one. Let us 
consider first the spectroscopic results. In the best case (all redshift bins over an area of 1000 
deg^ survey) the errors on wq and wi are 0.12 and 0.16, respectively. Marginalizing over G 
we obtain 0.19 and 0.23: the information in G reduces the errors by roughly 30% . The error 
on Qrn goes from 0.017 including G to 0.03 marginalizing over G. 

For the photometric case {6z = 0.02) the best case of Table IV gives a^o = 0.32 and 
cr^j = 0.39. A hybrid survey could reduce the errors to 0.23 and 0.23, respectively. Reducing 
the density of the z ^ 1 surveys leaves almost unaltered the constraints in the spectroscopic 
case (since the ratio nP/{l + nP) in Eq. (fT^ is already close to saturation) while weakens 
the photometric constraints by a few percent (see Table IV). 

The relative error on the Hubble parameter h = (um/^moY^'^ is dominated by the relative 
error on and is therefore approximated by an^o/^^mo- We find Ah/h 3 — 7%. 

It is interesting to comment also on the inclusion of SNIa data. The supernovae constrain 
Wo to a much higher precision than Wi. Multiplying the baryon oscillations contours with the 
SN limits one sees that the advantage of spectroscopy reduces considerably, at least below 
200 deg^, almost regardless of 5z, as can be seen in Fig. 7. 
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We compare our baryon oscillations forecasts also with the present SN data. In Fig. 8 
we show the constraints from the current set of supernovae of Riess et al. (2004; gold sam- 
ple), along with the constraints for a 200 deg^ survey (marginalizing over G) with reference 
cosmology centered on the Riess et al. (2004) best fit, Wq ~ — 1.3,Wi ~ 1.5 and Vtrno = 0.27. 
This plot shows that the size and orientation of the Fisher contour regions depends in a 
significant way on the underlying reference cosmology, as already remarked in SE. It also 
reveals that, in cases as this, the spectroscopic surveys maintain their advantage even using 
the SN information. 



Table 


III. Spectroscopic 


: surveys 


Area (deg^) 












z ■■ 


= 1 






100 


0.0053 


0.037 


0.39 


0.56 


200 


0.0053 


0.033 


0.34 


0.48 


1000 


0.0052 


0.024 


0.22 


0.31 




z ■■ 


= 3 






100 


0.0051 


0.045 


0.32 


0.36 


200 


0.0050 


0.043 


0.31 


0.33 


1000 


0.0041 


0.037 


0.25 


0.25 


all surveys 


100 


0.0051 


0.034 


0.26 


0.32 


200 


0.0049 


0.029 


0.21 


0.26 


1000 


0.0040 


0.017 


0.12 


0.16 
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Table IV. Photometric {6z = 


0.02) surveys 


Area (de^ 








z = 1 


100 


0.0054 0.046 


0.45 


0.72 


200 


0.0054 0.046 


0.42 


0.63 


1000 


0.0053 0.043 


0.35 


0.45 


z = 3 


100 


0.0053 0.046 


0.45 


0.70 


200 


0.0053 0.046 


0.42 


0.62 


1000 


0.0052 0.045 


0.35 


0.45 


all surveys 


100 


0.0053 0.045 


0.42 


0.62 


200 


0.0053 0.044 


0.39 


0.54 


1000 


0.0050 0.039 


0.32 


0.39 


hybrid surveys 


100 


0.0052 0.044 


0.31 


0.34 


200 


0.0050 0.042 


0.30 


0.31 


1000 


0.0041 0.033 


0.23 


0.23 




density/10 at 2; ^ 


1 




100 


0.0053 0.045 


0.44 


0.67 


200 


0.0053 0.045 


0.41 


0.59 


1000 


0.0051 0.040 


0.34 


0.42 
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Figure 2. Relative errors for H,D,G for the various surveys. The legend label x,y means a survey with redshift error x/100 
and area y ■ 100 deg'^. The dot at 2 = 0.3 is the SDSS survey. 



5 SCALAR FIELD MODEL 

The perfect fluid DE model of the previous section is a useful toy model but has no theoretical 
motivation and is incomplete in the sense that it is only a description of the background 
dynamics at low redshifts. In this section we perform the analysis by considering an example 
of a microphysical description of DE based on a scalar field model. 

We adopt to this purpose one of the first and simplest scalar field model proposed in 
literature, the inverse power-law potential (Ratra & Peebles 1988, Steinhardt et al. 1999) 
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Figure 3. Effect of the growth function for spectroscopic surveys: confidence regions marginahzing over G (dotted curves) and 
using the information in G (full curves). The plots are for the survey at z = 3 alone (upper panels) and for 2 = 1 alone (lower 
panels). In the color version of the figures lines referring to G-marginalization are blue, those referring the the inclusion of G 
are in red. 



V = . (25) 

The cosmological equations for this model are conveniently written in the e-folding time 
a = log a = — log(l + z) as (the prime denotes derivation with respect to a) 

0" + 0'(2 + ^) + aV^ = 0, (26) 
Id 

= H'^n^oa-' + ^i^<P"H' + V), (27) 

where = SttG, Hq^ = SOOOMpc, and ^k,p denote the kinetic and potential energy, 
respectively, of the scalar field. It is convenient to introduce the notation x'^ = = Q.k, 
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Photonietric surveys (5z = 0,02j 





Figure 4. Effect of the growth function for photometric surveys: confidence regions marginaUzing over G (dotted curves) and 
including G (full curves). The plots are for the survey at 2 = 3 alone (upper panels) and for z = 1 alone (lower panels). 



(28) 



y"^ = 1^ = fip and h = H/ Hq (which should not be confused with the dimensionless Hubble 
parameter). Eq. (jTFjl becomes then (Amendola 2000) 

_3 

^ = 1 2 2- 

Then the equations can be written as a dynamical system as 

X' = -^X(x2-y2_l)_^^2(l+n)/n^2/n^ 

y' = W+")/"/.2/« + ^y(l + x^-y^), 

h' = -|(3 + 3x^-31/2), 
where we introduced the dimensionless quantity 

-2 - [3 



(29) 



= 7^{—)~~A-^^ Ho 



n K , n + 2 

where 0o is the present value of the scalar field and Qpo its present potential density. The sys- 
tem can be solved numerically to derive x{z),y{z) for each value of the parameters n, flmo, 
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Martina lization on G 




Figure 5. Confidence ellipsis for the combined data CMB,SDSS, surveys at z 1, 3, for three redshift errors, 0, 0.02 and 0.04, 
inside to outside. Here and in all subsequent plots the redshift errors are indicated as 100 ■ Sz. The central dot is the reference 
cosmology. The bottom right panel refers to the hybrid survey discussed in the text (spectroscopic plus photometric): here 
the contours are for surveys of 1000, 200 and 100 deg^, inside to outside. The shaded areas are the confidence regions for 400 
simulated SNIa distributed according to the reference cosmology up to 2 = 1.5, at 68%, 95% and 99% c.l., inside to outside, 
marginalizing over Hq (uniform prior) and over ilrnO with a Gaussian centered on the reference value and variance 0.1. The 
dotted contour (at 68% c.l.) assumes instead a 10% Gaussian error on Hq for the SNIa. 



but this is computationally cumbersome and hides interesting features of the solutions. In 
Steinhardt et al. (1999) it has been shown in fact that such a model possesses a tracking 
solution, defined as a trajectory on which many initial conditions converge, as shown in Fig. 
9. Once a solution reaches this trajectory its properties are determined solely by the slope n 
and the present value flmo- As long as x, ?/ -C 1 the tracking trajectory has ~ — 2/(2 + n); 
when later on DE takes over, y tends to unity and — > —1. Since the variation of is 
relatively small, to a first approximation the inverse power-law model for n <^ 1 resembles a 
perfect fluid with ~ — 2/(2 + n) and Wi ~ 0. 

In the Appendix we derive a semi-analytical fast approximation to the tracking trajectory 
that interpolates correctly between the two limits of w^. This approximation will be used in 
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Figure 6. As in Fig. 5 but including the dependence of the growth factor G on the DE parameters. 



what follows; it is found to be precise to better than 1% on VL^n^z) and to better than 5% on 
w^{z) across all the relevant range. Further, we find that for the growth function the same 
approximation d\ogG{z) / da = Q'^la) adopted for ACDM remains valid, with deviation 
from the numerical result less than 0.1% in the range z G (0, 3) putting 7 = 0.575 for n ~ 
and 7 = 0.625 for n ^ 0.5. For all practical purposes the dependence of 7 on n can be 
neglected and we put 7 = 0.6 in the whole interesting range. 

We project now the Fisher matrix in the subspace 
dn = i^moh^, ^mo, n) . (30) 

The results are shown in Fig. 10 (reference cosmology n = 0.05, practically indistin- 
guishable from ACDM ) projected on the plane Qrno,n, and in Table V. As expected, the 
SN contours are elongated towards the contours of HD (see the contour plots in Fig. 11). 
The contours of the Fisher matrix are instead slanted towards the direction of the iso-G 
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Wo 



Figure 7. Contours at 68% of the combined likelihood of supcrnovae and surveys of 200 deg^ with redshift errors 5z = 
0,0.02,0.04, inside to outside (including G). 



lines: this shows that it is the information in G that drives the results. This optimizes the 
complementarity between SNIa and LSS. 

The tracking scalar field model has been tested towards SNIa and CMB. We analysed 
the recent Riess et al. (2004) SNIa data with the inverse power-law potential and found 
the constraint n < 0.9 at 68% c.l. [n < 0.5 assuming Q^o = 0.3 ± 0.05). On the other 
hand WMAP data allow for n < 0.94 (or n < 0.74 including the HST prior on the Hubble 
constant), at 68% c.l (Amendola & Quercellini 2003). From Fig. 10 and Table V we see that 
the baryon oscillations for a scalar field model can constrain Qmo to within 0.01 to 0.03, and 
the slope n to n <0.15 and n <0.53, depending on the survey strategy. In all cases this is 
quite better than the current and near-term SN and CMB constraint. Once again, most of 
the advantage of spectroscopic measures vanishes if one takes into account the SNIa contours 
(see Fig. 12) 

Here again we perform the comparison with the Riess et al. (2004) gold sample of su- 
pernovae (Fig. 13). It can be seen that the present SN data already cut most of the vertical 
elongation of the contours for the photometric surveys, leading to constraints on n which are 
similar to those from spectroscopic surveys. The spectroscopic errors on Qmo remain however 
competitive. 
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-2 



Figure 8. Here we compare the confidence regions of 200deg^ surveys witli redsiiift errors of 0,0.02 and 0.04 (marginalized 
over G) to the recent SN data of Riess et al. (2004) assuming fimO = 0.27 it 0.04 (Gaussian errors). For this plot the reference 
cosmology in the Fisher matrix coincides with the best fit to the SN {wq = — 1.3,uii = —1.5) . The SN contours are 68,90 and 
95% . The star marks the ACDM model. 



-0.2 



^ -0.6 




Figure 9. Equation of state of the scalar field model for several different initial conditions. The plot shows the tracking behavior: 



all curves converge to a trajectory which, at early times, has constant equation of state 



to to J, 



-1. Here n = 0.5 and Q„ 



: 0.35. 



2 + n 



and asymptotically converge 
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6 CONCLUSIONS 

Future deep and wide redshift surveys will have the chance to explore the little known 
landscape at z ~ 1 — 3 between the recent and early universe. This is the epoch at which 
dark energy begins to impel its thrust to the expansion: it is therefore of extreme interest 
for cosmology in order to gain hold of reliable data on its evolution. As it has been shown 
by the recent releases of high redshift SN data (Tonry et al. 2003, Riess et al. 2004) the use 
of standard candles will have much to improve before it will be able to set stringent limits 
on the dark energy equation of state evolution. Moreover, there are still no solid proposals 
for standard candles visible at redshifts significantly higher than z ~ 1.5. Other probes of 
the expansion djTiamics at larger distances (gamma-ray bursts, lensing) might be promising 
but depend on unknown physics and modeling. 

The use of the baryon oscillations in the matter power spectrum is on the contrary 
based on well-known phenomena that have already seen a spectacular validation in CMB 
observations. Moreover, it can provide a tomography of the universe expansion from z = 
down to z ^ 3.5, virtually without loss of precision. An intrinsic advantage of this method 
is the combined use of H, D and G to reduce the degeneracies of each of these quantity. 
We have shown that the inclusion of G reduce the errors by 30% roughly. We found that 
a 200^ survey with absolute redshift error of Sz = 0.02 can limit Wq,Wi by 0.39 and 0.54, 
respectively, assuming wq = —l,wi = as reference values. A spectrocopic survey can reduce 
the uncertainties to 0.21 and 0.26. If the dark energy is modeled as a scalar field with inverse 
power-law potential n, the limits are n < 0.26 (spectroscopic redshift) and n < 0.40 (redshift 
error 6z = 0.02), along with an estimation of Qmo to better than 3%. Inclusion of supernovae 
data further reduce the errorbars due to the different direction of degeneracy. The reduction 
is found to be particularly efficient in the case of photometric redshift measurements. 
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APPENDIX. THE FAST-TRACK APPROXIMATION 

A tracking solution, as any other trajectory, is defined in a 3D phase space by two relations of 
the form ei_2 = fi,2{x, y, h) where ei,2 are constants. There must exist therefore two constants 
of motion of the system (jSHI)- In the original derivation of the tracking solution (Steinhardt 
et al. 1999) it was found that the equation of state 

x'^ + y^ 

is almost constant during the tracking. The second constant of motion is 

q = y{2+-)/"/,2/n _ (32) 

Then, by derivating, we obtain for the tracking solution 



qt = - — Vl - w 



2 

2/i ' ' " ' 



2 + n{y^-x^) 

wt = — • 33 

n + 2 

The second relation says that w, q are constant only for y"^ — x"^ <^ 1: this is the regime in 
which the tracking behavior was originally found. In this case 

w = — . (34) 

2 + n ^ ' 

However, since near the present epoch x^,?/^ are no longer negligible, it is necessary to 
improve the accuracy over the solution (jH^ . Substituting w from Eq. in we obtain 
(assuming w > —1) 
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Figure 10. As in Fig. 5 but for tlie scalar field model. Here the SNta confidence regions are not marginalized over QmO- 



w{y,n) 



l-2y^- J9- 24(1 + n)/n + 16(1 + n)2/n2 - 4y'^ + %4 



A/n + 2 



From (jH^ we obtain 

y = (qh^/^)2T^ . 

Employing (|28p we obtain finally 



(35) 



(36) 



2/n 



- — ^Jl-w{y,ny 



l-yH , 1 J 

^ l—'w{y,n) I 



n+2 



(37) 
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Figure 11. The plot displays the contours of D,HD,G at 2: = 1 (left column) and 2 = 3 (right column). Each function is 
normalized by the value at the reference cosmology, marked by the white dot. The contours are for steps of 0.05, increasing 
from black to white 



© 0000 RAS, MNRAS 000,0-?? 



Constraints on perfect fluid and scalar field dark energy models from future redshift surveys 25 



0.8 



a 0.6 



0.4 



0.2 



400 SN+200 deg^ survey 
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Figure 12. Scalar field; contours at 68% of the combined likelihood of supernovae and surveys of 200 deg^ with redshift errors 
Sz = 0,0.02,0.04, inside to outside (including G). 



The parameter fj, has to be fixed so that at the present + = 1 — flmo ■ This gives a 
relation = /i(n, ^mo)- We found by iterative integrations that a very good fit to is 



, „ _ 0.00036 0.0033 + 0.0088(]^o 
exp[— 0.53 1 



0.13 - 1.16^]„o + 0.41fi^o - l.mVLl 



n 



+ 0.31n + 2.43fi 



mO 



With this substitution, Eq. (jHTjl becomes an imphcit algebraic relation that, solved numeri- 
cally for y, gives |/(a; n, VLmo) and all related quantities (a;(t), w{t), h(t)). 
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Table V. Results for scalar field. 


Area (dej 


I ) '^UJm 






5z = {) 


100 


0.0032 


0.018 


0.32 


200 


0.030 


0.016 


0.26 


1000 


0.0022 


0.010 


0.15 


5z = 0.02 


100 


0.0034 


0.030 


0.44 


200 


0.0033 


0.028 


0.40 


1000 


0.0031 


0.021 


0.30 


6z = 0.04 


100 


0.0034 


0.032 


0.47 


200 


0.0034 


0.030 


0.43 


1000 


0.0033 


0.026 


0.34 


hybrid surveys 


100 


0.0032 


0.023 


0.39 


200 


0.0031 


0.020 


0.34 


1000 


0.0026 


0.014 


0.23 
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Figure 13. Here we compare the confidence ellipsis of 200deg^ surveys with absolute redshift errors of 5z = 0,0.02 and 
0.04(inside to outside), with the confidence regions for the scalar field model with the SN data of Riess et al. (2004) For this 
plot the reference cosmology in the Fisher matrix coincides with the best fit to the SN [QrnO = 0.3, n = 0). The SN contours 
are 68,90 and 95% . 
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